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Abstract 

Using piezoresponse force microscopy on epitaxial ferroelectric thin films, we have measured 
the evolution of domain wall roughening as a result of heat-quench cycles up to 735°, with the 
effective roughness exponent £ changing from 0.25 to 0.5. We discuss two possible mechanisms for 
the observed £ increase: a quench from a thermal 1-dimensional configuration, and from a locally- 
equilibrated pinned configuration with a crossover from a 2- to 1-dimensional regime. We find 
that the post-quench spatial structure of the metastable states, qualitatively consistent with the 
existence of a growing dynamical length scale whose ultra slow evolution is primarily controlled by 
the defect configuration and heating process parameters, makes the second scenario more plausible. 
This interpretation suggests that pinning is relevant in a wide range of temperatures, and in partic- 
ular, that purely thermal domain wall configurations might not be observable in this glassy system. 
We also demonstrate the crucial effects of oxygen vacancies in stabilizing domain structures. 

PACS numbers: 77.80.Dj, 68.35. Ct,77.80.Fm 
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I. INTRODUCTION 



Domain walls separating differently polarized regions in ferroelectric thin films present 
a powerful model system in which the characteristic roughening and complex dynamics 
resulting from competing elastic and pinning forces can be readily accessecR Under- 
standing this behavior is key to the physics of many diverse disordered elastic systems 
(DES)P^ and of significant technological interest for memory and electromechanical fer- 
roelectric application^^, and devices using the domain walls as nanoscale functional 
component ^ 13 1 14 1 , in which control over the stability and growth of domain structures is 
of paramount importance. 

Theoretically, although equilibrium DES properties are relatively well understood, much 
less is know about the out-of-equilibrium behavior of these systems. A major challenge is 
understanding the non-steady slow dynamics associated with agingpHTD A specially inter- 
esting realization of such out-of-equilibrium phenomena is provided by quenches, in which a 
parameter is abruptly varied. The study of such quenches is a challenge, not only for clas- 
sical systems but also for quantum system^El. p Q r classical systems, studies of interfaces 
subjected to quenche^^J h ave shown aging of the interface and a long term memory of the 
initial configuration. 

On the experimental side, studies of ferroelectric domain walls in the DES framework have 
so far focused mostly on equilibrium properties, with switching current and dielectric mea- 
surements showing domain wall pinning in bulk ferroelectrics 25 , and random-field disorder 
inferred for relaxors 26 , while piezoforce microscopy (PFM) studies of domain wall roughen- 
ing and creep reported strong dipolar interactions and random bond disorder in epitaxial 
ferroelectric thin film^HSU _ u n iik e their theoretical counterparts, these systems also show 
pinning by individual strong defect^, and local variations of both disorder strength and uni- 
versality clasd^, further complicating the approach to non-equilibrium behavior. Thermal 
effects in bulkpl and thin film material^ include increased roughness and higher domain wall 
mobility upon heating, and transition to domain wall depinning at lower temperature^. 
However, a detailed nanoscale analysis of the complex behavior associated with a thermal 
quench of these systems remains an open question. 

In this paper we report on PFM measurements of domain walls in ferroelectric 
Pb(Zr .2Ti .8)O3 (PZT) thin films, heat-quench cycled up to 735°C, in which the value 
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of the roughness exponent £, reflecting the system dimensionality and the nature of the 
disorder, is extracted from an analysis of the domain wall position, and evolves from ^0.25 
to ^0.5. Despite the similarity of the measured exponent with the thermal = 1/2 in one 
dimension, we show that the observed behavior cannot result from the quench of an initial 
high temperature thermal configuration, as the observed roughness values are much higher 
than would be expected in this scenario. Rather, it is clear that disorder pinning continues 
to dominate the behavior during the quench, with one possibility being that the domain 
walls cross from a 2- to 1-dimensional regime. We also follow the complex thermal evolu- 
tion of circular nanoscale ferroelectric domains, in which the combined effects of disorder 
pinning, line tension and a preferential polarization orientation can be observed, and the 
crucial effects of oxygen vacancies in stabilizing domain structures are highlighted. 

II. MATERIALS AND METHODS 

Thermal quench studies were carried out on c-axis oriented PZT films, epitaxially grown 
on single crystal (001) Nb:SrTi03 substrates at ^500° C by off- axis radio- frequency mag- 
netron sputtering, in an Ar-02 process gas mixture^, with high crystal and surface quality 
confirmed by x-ray diffraction and atomic force microscopy (AFM). In these films, the po- 
larization vector along the c-axis can be locally switched by a biased AFM tip, using the 
substrate as a ground electrode, and the resulting ferroelectric domains imaged by PFlVp^. 
To study the effects of heat-quench cycling on domain wall roughness we wrote 10 x 10 
/xm 2 arrays of linear 1 fim wide domains with ±12 V alternately applied to a scanning 
AFM tip. To explore the interaction between heating, line tension and disorder pinning, 
we wrote arrays of nanoscale circular domains with ±12 V pulses of different duration in a 
uniform, oppositely polarized region. The films were heated in air to progressively higher 
temperatures for fixed time intervals, then quenched to room temperature on a large copper 
block. Ambient temperature PFM images of the evolving domain structures were acquired 
after each heating-quench cycle. In addition, to investigate the role of oxygen vacancies on 
domain stability, films grown under the same conditions but cooled in either process gas, 
Ar, or O2 were compared. 

As in our previous workpS to extract the correlation function of relative displacements 
B(L) for a given PFM phase image of linear domains we first isolated a single domain 
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wall, then binarized the image using a cutoff extracted from the midway point between 
the phase contrast values for the two polarization orientations. The elastically optimal flat 
domain wall position was determined from a least-squares linear fit to the binarized image, 
which finally allowed the extraction of the relative displacements, and their correlation at 
different length scales. To compensate for possible effects of scanner drift, for each set 
of walls PFM measurements taken in an "up" vs "down" scan direction were averaged. 
Although ambient-condition walls after writing were generally well defined, repeated heat- 
quench cycles, especially close to Xc, led to the appearance of multivalued features such as 
bubbles or overhangs. In some cases, such features could also be present even immediately 
after writing at room temperature. Bubbles, being very small, well defined domains of 
opposite contrast near but not coincident with the position of the primary domain wall, 
could be removed from the binarized image by hand. For overhangs, the algorithm used to 
calculate relative displacements would simply result in a renormalized single value "effective 
position" of the domain wall. Therefore, when such features became prevalent, or when the 
films switched to the preferred polarization direction upon heating near Xc, quantitative 
analysis could no longer be carried out, as indicated by the symbol "x" in Table. [T| 



III. DOMAIN WALL ROUGHNESS 



Generally, an elastic manifold in the absence of disorder and thermal fluctuations (OK 
temperature) would take an ideally flat configuration in order to optimize its energy. Ran- 
dom variations in the potential landscape, whether due to disorder or thermal fluctuations, 
allow further optimization as the manifold wanders between particularly favorable regions 
in the potential, giving rise to a characteristic roughening. In the case of single manifold, an 
interface such as a domain wall, contact line, or fracture, fully in equilibrium with its random 
potential landscape, the resulting geometrical fluctuations from an elastically optimal flat 
configuration are present at all length scales, with the system showing mono-affine scaling 
properties (for a more detailed presentation, see and references therein). The key quantity 
describing this roughening is the correlation function: 

B(L) = {[u(z + L)-u(z)}i)<xL 2 < , (1) 



4 



known as the roughness function. Here u(z) are the transverse displacements along the 
longitudinal coordinate z with respect to a flat configuration, measured a distance L apart 
on the interface, and (...) and ^ the thermal and disorder averages, respectively. The value 
of the roughness exponent ( depends on the dimensionality of the system, the nature of 
the disorder potential, and the range of the elastic interactions. Extracting ( not only 
yields this important information, which may be used to determine the universality class 
of the disorder and the dominant pinning defects present in the system, but also allows 
a full scaling prediction of the complex non-linear creep response of such systems to small 
driving forced. We have previously shown for linear domains created artificially by straight- 
line scanning with a biased AFM tip that, as a result of this writing process, the domain 
walls are essentially flat at large length scales, with roughening observed only at short 
length scales, probably where the stray fields during the writing process allow some small 
accommodation to the potential landscape^. The experimentally obtained B(L) thus shows 
a power-law growth at short length scales, from which the roughness exponent ( can be 
extracted, followed by a saturation at large length length scales. 

As previously reported irP2 linear domains in process-gas-cooled films show very high 
stability upon heating, with increased roughness only at small length scales (Fig. [l](a-f)). 
At 625~785°C (depending on the sample), close to the Curie temperature Tc of the films 
determined by x-ray diffraction^, extensive polarization switching occurs and the writ- 
ten domains eventually disappear as the region reverts to its as-grown monodomain state. 
Quasi-identical behaviour was observed in the Ar-cooled samples (Fig. [T](m-o)). In the 
02-cooled films however, the domain structures showed much lower thermal stability. In 
these films, expected to show decreased oxygen vacancy densities compared to process-gas 
or Ar-cooled samples, increased roughening was found even for short (20 min.) heating inter- 
vals, with domain disappearance observed already at 250-350°C (Fig. [T](j — 1) ) . Nonetheless, 
x-ray measurements confirmed sample Tc in the 600-700°C range, similar to those of the 
process-gas-cooled sampled. 

To quantitatively examine the roughness evolution, for each sample we extracted the 
average roughness function B(L). For the relatively flat, as- written domain walls before 
thermal cycling, as in our previous studied, we find a power-law growth of B(L) at short 
length scales followed by saturation at L* ^25-100 nm (comparable to the ^40-100 nm 
film thickness), with £>(L*) ^100-1000 nm 2 in the different films. After heating, the power- 
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FIG. 1: PFM images of linear ferroelectric domains in a process-gas-cooled 50 nm PZT film (a) 
immediately after writing and (b-g) after thermal cycling to the indicated temperatures (°C). 
Similar measurements (j-1) in an 02-cooled 100 nm film and (m-o) an Ar-cooled 40 nm film, where 
higher oxygen vacancy densities appear to stabilize the domain walls. Vertical 1 /im scale bars 
are shown in all images, (h) Average B(L) = ([u(z + L) — u{z)] 2 ) extracted after thermal cycling 
to the indicated temperatures (°C) for all the domain walls in the 50 nm film, (i) The power- 
law growth region of B{L) for same measurement (identical color code for the temperatures). The 
extrapolated fit to the power-law growth region of B(L) (first 20 points) with ln[£> (L)] oc 0.890 ln[L] 
is shown in red, and the fit to the saturation region of B(L) in aquamarine for the 685° data set. 
The intersection of these fits allows the saturation length scale L* to be extracted. The dashed 
vertical line in (h) and (i) indicates the length scale L equal to the film thickness. 

law- growth region of B{L) extends to L* ^300-700 nm and the saturation value increases 
(£>(L*) ^1000-10000 nm 2 ), consistent with the observed increase in domain wall roughness 
(Fig. [j](h,i)). From the power-law growth region of B{L) (Fig. [jjl)), we extract a value 
for the roughness exponent £, characterizing the static equilibrium configuration of the 
domain wall in the random manifold regime, where B(L) oc L 2 ^. As shown in Table, [i] 
for domain walls in different films, heating increases the value of ( from ^0.25 at ambient 
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conditions and lower temperature thermal cycling to ^0.5 at higher temperatures and near 
domain disappearance^. In 02-cooled films, higher roughness and increased ( values are 
seen already at relatively low temperatures. 

TABLE I: Thermal evolution of the roughness exponent ( for PZT films cooled in process gas 
(PG), argon (Ar), and oxygen (O2) after long (A: 1-3 hrs., P.P.) and short (B: 20-30 min., X.H.) 
heat-quench cycles. - data set missing or noisy, x temperatures beyond domain disappearance, 
multivalued wall roughness configuration. * 225° C measurement. 



A) T(°C) 


25 


110 


200 


285 


385 


485 


585 


685 


50 nm PG 


0.33 


0.43 


0.42 


0.48 


0.54 


0.50 


0.51 


0.45 


91 nm PG 


0.23 


0.22 


0.35 




0.37 


0.39 


X 


X 


52 nm PG 


0.26 


0.26 


0.18 


0.23 


0.26 


0.31 


0.45 


X 


B) T(°C) 


25 


200 


300 


400 


475 


500 


575 


600 


50 nm PG 


0.10 


0.15 


0.31 


0.30 


0.27 


0.42 


X 


X 


40 nm PG 


0.24 




0.29 


0.29 




0.24 


0.42 


0.50 


40 nm Ar 


0.27 


0.18 


0.17 


0.27 




0.22 


0.38 


0.41 


100 nm 2 


0.15 


0.45* 


0.47 


X 


X 


X 


X 


X 



IV. DISCUSSION 

To understand these data, we first consider the thermal roughening of the linear domain 
walls within the theoretical framework of DES. Our previous studies have shown that at 
ambient conditions with no applied electric field, thermal excitations alone are not enough 
to overcome the energy barriers between different metastable states in the duration of the 
experiment, although the electric field applied during writing allows domain walls to accom- 
modate at short length scales to the surrounding random potential landscape. A power-law 
growth of B{L) is thus observed only at length scales smaller than the film thickness, where 
the domain wall behavior agrees well with theoretical predictions for 2-dimensional inter- 
faces pinned by weak, collective random bond disorder and in the presence long-range dipolar 
£ orce ^28|37j ! Given the finite thickness of the sample, much smaller than its lateral dimen- 
sions, it is reasonable to expect that in a weak collective pinning scenario, were the walls to 
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fully equilibrate, 1-dimensional behavior would be evident at large length scales. Although 
it is difficult to estimate the exact length scale L x for such a 2-dimensional to 1-dimensional 
crossover in a disordered system, we can expect it to be of the order of the sample thick- 
ness. This would be exactly the case for purely thermal roughening. The increased energy 
provided to the samples on thermal cycling to progressively higher temperatures is one way 
to promote further equilibration of the domain walls, allowing them to accommodate to the 
surrounding potential landscape at length scales greater than the film thickness, where the 
1-dimensional nature of the domain walls should become evident. 

If thermal roughening dominates in the high temperature 1-dimensional regime, the walls 
should be characterized by a roughness exponent = 1/2, and slowly age to the equilibrium 
( eq configuration when quenched. For an increasing time t after the quench, a growing length 
Ldyn(t) should thus separate the equilibrated short length scales (with exponent ( eq ) from 
the large scales conserving a memory of the roughness at the quench (with exponent (^), 
analogous to a quench from an initially flat configuration^. We confirmed such a scenario 
in Langevin dynamics simulations of a one-dimensional interface in a random medium as 
detailed in Refl2Q[ but this time quenching an initially equilibrated high-temperature config- 
uration. The results indeed show the quasi-freezing of this initial configuration at the largest 
length scales, with thermally activated dynamic evolution distinguishable only at very small 
length scales (Fig. [2](a)). This behaviour is also reflected in the evolution of £>(L, t) after the 



quench (Fig. [2](b) ): at large lengths L > Ld yn (t) we observe B(L,t) ~ L 2C,th , while at small 

lengths the line tends to approach the equilibrium roughness exponent ( eq (with ( eq = 2/3 
for an elastic line with short range elasticity (random bond disorder)). 

If the barriers are large, as the high stability of ferroelectric domain walls indicates, the 
length Ldyn(t) can grow extremely slowly, and most of the observable interface could well be 
characterized by the exponent = 0.5. This scenario seems consistent with the exponent 
( ~ 0.5 observed at length scales up to 10 times greater than the thickness of the film. 
However, if we consider an elastic 1-dimensional interface with energy H = | / dz(Vu) 2 , we 
obtain its thermal roughening at temperature T as 



where the elastic constant c in the case of the domain wall is the product of its energy per unit 
area (Jdw an d the film thickness df and ks is the Boltzmann constant. Taking (Jdw — 132 
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FIG. 2: Numerical simulation of an overdamped elastic string in a random medium undergoing a 
temperature quench T\ — >> T2 < T\ at t = 0. (a) From bottom to top: Typical evolution of the 
elastic string each 10000 steps. The bold line is the configuration at the quench, (b) From top 
to bottom curve: Evolution of the disorder- averaged correlation function of relative displacements 
B(L,t) after the quench (symbols). The roughness exponent Qh — 0.5 describes the geometry of 
the string at large length scales (solid line). 

mJm" 2 from ab-initio calculations for PbTiOj^, we can therefore calculate the gradient 
S = ^2c~ = 2aowd f 7 an< ^ ^ e rou g nness B(L*)th expected for the thermal roughening scenario 
at the length scale corresponding to the experimentally observed saturation length L*. These 
values can be directly compared to the experimentally obtained values for the roughness 
5(L*) exp and the gradient S exp = B{L*) exp /L* in the power-law-growth region where B{L) oc 



L 2C> . As shown in Table. II for the 50 nm PZT film of Fig. 1, the experimental B{L) values 
are many orders of magnitude higher than those expected for a purely thermal roughening 
of a 1-dimensional interface. The observed behavior at large length scales therefore cannot 
be simply a quench of the high temperature configuration. Rather, pinning by disorder in 
the film continues to be the dominant effect, even during the quench. 

One possible scenario is that of a dimensional crossover dominated by disorder: in equi- 
librium with a random bond disorder potential, below a typical length scale L x , compa- 
rable to the film thickness, the domain walls should act as 2-dimensional elastic sheets 
with ( 2 d — 0.2-0.3, and as 1-dimensional elastic strings with d D = 2/3 for higher length 
scales (see Figj3](d)). If this crossover is smeared out, it is not improbable that an effective 
Cx : (2D < Cx < (id is observed experimentally. According to the theory of DES 20 , the 
locally equilibrated regime just after the quench should extend up to a length-scale that 
grows with the heating-time t and the heating temperature T, as the dynamical length 
Ld yn (t) ~ [(T/Uq) log(t/t )] 1 ^ 6> 7 with 6 a positive exponent, Uq a characteristic temperature, 
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TABLE II: Comparison of the roughness function B{L) and its gradient S — ^ff- for a thermally 



roughened 1-dimensional interface with the experimentally obtained gradient of B(L) oc L 2 ^ in a 
50 nm PZT film. 



T(K) 


383 


473 


558 


658 


758 


858 


958 


5 x 10" 4 (nm) 


4.0 


4.9 


5.8 


6.9 


7.9 


9 


10 


L* (nm) 


58 


121 


72 


180 


199 


200 


232 


B(L*) th (nm 2 ) 


0.023 


0.060 


0.042 


0.124 


0.158 


0.179 


0.232 


B(L*) exp (nm 2 ) 


275 


1593 


1143 


2805 


3065 


4089 


4343 


S exp (nm) 


4.8 


13.1 


15.8 


15.5 


15.4 


20.5 


188 



and t a characteristic microscopic time. We then expect to observe equilibrium power 
law behavior (possibly including the equilibrium dimensional crossover) below Ld yn (t) and 
a memory of the flat initial condition, imposed by the linearly scanning AFM tip during 
writing, above it. The theory thus predicts a temperature dependence of the effective rough- 
ness exponent, from ( x ~ ( 2 d at low temperatures, T : Ld yn (t;T) ~ L x (see Figj3](b)), to 
Cx ~ (id at higher temperatures T : Ld yn (t]T) 3> L x (see Figj3](c)) where the power law 
behavior is effectively dominated by the equilibrium 1-dimensional regime. The results dis- 
played Fig(l](e,l) qualitatively agree with this theoretically expected scenario, and in partic- 
ular, the results in tableware consistent with the expected heating-temperature dependence 
of the effective roughness exponent. The theory also predicts, at a fixed temperature, a 
slow (logarithmic) time dependence of Ld yn (t) that should translate into a heating-time de- 
pendence of Cx, from ( 2 d at short heating times t : Ld yn (t]T) ~ L x , to (id at large times, 
t : Ld yn (t) 3> L x , as summarized in Figj3] . In this respect, the results shown in Table [T] for 
the 50 nm PG-cooled samples in (A) and (B) seem to be also qualitatively consistent with 
the theory. 

Moreover, the S exp values in table II show no clear linear temperature dependence, a 
fact that is also qualitatively consistent with large length scales dominated by disorder, and 
which we can further analyze through the following estimates. Phenomenologically, in the 
case of disorder-dominated roughening of an elastic interfaces] ? we can consider in the 
2-dimensional regime at short length scales that B{L) = r 2 {L/ Lc) 2C,<k2D \ where Tj is the 
correlation length of the disorder potential and Lq the Larkin length (minimum pinning 



10 



length, below which the interface would behave purely elastically) . In the 1-dimensional 
regime at large length scales, this relation would become B(L) = r 2 char {Lj L X ) 2 ^ 1D \ where 
the characteristic length in this case would be the crossover length scale L x and r c h ar = 
r f (L x /L c ) c{2D) . Taking L c ~ 0.2 nm from|29l ((2D) = 0.25, ((ID) = 0.5, and L x - d f 
for the 50 and 91 nm samples, we find that the experimentally obtained B(L) values would 
correspond to r$ ~ 4 nm in the 2-dimensional regime (ambient conditions) and r$ ~ 4-8 nm 
in the 1-dimensional regime (thermally cycled). For reference, the 10 18 -10 19 cm -3 impurity 
densities extracted from current-voltage measurements in PZT filmd^ correspond to an 
inter-defect spacing of ~ 4.5-10 nm, and would thus qualitatively agree with the extracted 
disorder correlation lengths. 

An additional important consideration, especially at increasing temperatures, is the exact 
nature of the disorder. Although surface deterioration is seen upon heating, especially at 
the highest temperatures, x-ray characterization of the samples demonstrates continued 
high crystalline quality 33 . However, the mobility of point defects, such as oxygen vacancies, 
increases strongly with temperature. From the comparison of PG-, Ar- and 02-cooled films, 
it appears that the primary defects responsible for the domain wall pinning are indeed oxygen 
vacancies. These defects can appreciably lower their energy by associating with a domain 
walP^, so the probability of a migration and segregation of defects at the domain walls upon 
thermal cycling, effectively increasing the strength of the pinning potential and stabilizing 
the domain structures against large-scale thermal roughening, is very high. In this respect 
it is worth noting that in the theory of DES an analogous "slow structural relaxation" 
mechanism was very recently shown 42 to provide a minimal model for predicting pinning 
enhancement or aging of the static friction force, and the phenomenon of "aftershocks" or 
static-relaxation-induced instabilities in the interface. This is interesting as these phenomena 
are both absent in the standard DES models with truly quenched disorder, but they are for 
instance observed in solid friction experiments. 

Another possibility is that the thermal cycling steps, by removing the water layer normally 
present on the film surface in air at ambient temperature, significantly affect the electrostatic 
boundary conditions, and thus the effective potential landscape experienced by the domain 
walls. In this case, the higher roughness and increased ( values could reflect the configuration 
adapted by the walls under these new conditions, and frozen-in after a quench from high 
temperature, possibly in parallel to a dimensional crossover. Further investigation would be 
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FIG. 3: Disorder dominated scenario as a function of length scale, (a) At ambient temperature 
before heating, the as-written, initially flat domain walls are only able to locally equilibrate at 
short length-scales, (b) and (c) Heating facilitates thermally activated glassy relaxation and longer 
length scales become locally equilibrated as the heating temperature (or heating time) increases, 
(d) If fully in equilibrium (reference state), the system displays a 2-dimensional to 1-dimensional 
crossover at length scale L x , with no memory of the written state. This may be reflected by 
the effective roughness exponent ( e ff in the correlation function B(L) ~ L 2C,e ff of the quenched 
metastable states, as in (b) and (c), evolving from the 2-dimensional value (id, when Ld yn < L x , 
towards the 1-dimensional value (id when L^ yn ^> L x , especially if the crossover length scale 
is smeared out. In the Id regime the "dynamical length" L^ n (t;T) is expected^ to grow as 
Ldyn(t;T) rsj [(T/Uq) logft/to)] 1 / ', with 9 a positive exponent, Uq a characteristic temperature, and 
to a characteristic microscopic time. This prediction agrees qualitatively with the experimentally 
observed increase in £ and L* values with heating time and temperature. 

needed to asses the potential importance of this effect and also to quantitatively test the 
scenario described above, summarized in Figj3) 



V. DOMAIN SWITCHING AND GROWTH 

The crucial role of disorder in stabilizing the domain walls in a given configuration is 
further highlighted by the behavior of nanoscale circular domains followed in parallel mea- 
surements during the same heat-quench cycles. In these measurements, we note the strong 
polarization switching asymmetry: Pdown domains written with positive tip voltage in a 
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FIG. 4: PFM images of nanoscale ferroelectric domains written at ambient conditions in the 91 
nm PZT film (P UP as grown) with (a) 1 s, + 12 V, (c) 10 s, -10 V, (e) 1 ms, +12 V pulses, and 
in a 52 nm PZT film (Pdown as grown) with (g) 100 ms, +12 V pulses. Domain evolution after 
thermal cycling (f,h) to 200°C and (b,d,i) 585°C. Vertical bar = 250 nm. All images are shown to 
the same scale. The average domain size as a function of the writing time at ambient conditions, 
and after thermal cycling up to the indicated temperatures (j) for the 91 nm film and (k) for the 
52 nm film. 

negatively pre-polarized area grow significantly larger for a given pulse duration (Fig. [3](a- 
f)), even in films Pup polarized as-grown. We relate this feature to the highly asymmetric 
device configuratiorJ^H. In PG- and Ar-cooled films at ambient temperature, both tip volt- 
age polarities give uniform, circular domains, which become irregular upon heating, with 
apparent bowing of the domain walls. In P UP as- grown films, Pdown domains decrease in 
size after thermal cycling, while in Pdown as-grown films a slight increase in domain size 
was observed. Some of the smaller domains (20-50 nm radius) of both polarities collapse 
even after low-temperature heating intervals (Fig. [I](e-f)). In Pdown as-grown films, signif- 
icant displacements and an "amoeboid-like" growth are observed for Pdown domains near 
Tc, with the structures eventually coalescing (Fig. [4](g-i)). In the 02-cooled films, even at 
ambient conditions Pup domains were unstable and Pdown domains were stable only for 
domain radii greater <~ lOOnm. 

This complex thermal evolution can be understood by considering the forces on the 
domain wall due to its curvature and the polarization orientation as well as disorder and 
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thermal effects. A preferred monodomain state indicates an asymmetry of the ferroelectric 
double well potential, which promotes growth or collapse of domains depending on their 
polarization orientation. In addition, the "line tension" or elastic energy cost of circular 
domain walls promotes domain collapse, its effects diminishing with increasing domain size. 
Finally, disorder allows domain wall pinning, and therefore deformations from the elastically 
optimal circular geometry. The relative strength of these contributions during thermal 
cycling can be gauged from the degree of domain stability. Although polarization orientation 
and device asymmetry^ determine the domain size for a given writing time, for domains of 
comparable size, line tension determines stability. Larger P UP and Pdown domains remain, 
while 20^50 nm radius domains readily collapse on heating. Domain wall roughening is 
observed for all domain sizes, although the extent of the deformations increases with domain 
size. Disorder appears to be crucial in stabilizing the domain structures, in particular the 
smallest domains, and protects them against thermal effects. 

VI. CONCLUSION 

In conclusion, our studies show that in epitaxial ferroelectric thin films, domain wall 
roughening is governed primarily by disorder pinning, even upon heating to progressively 
higher temperatures, followed by a quench. However, the increased energy provided by the 
thermal cycling does allow accommodation of the artificially straight AFM- written domain 
walls to the surrounding disorder potential at progressively larger length scales, with the 
roughness exponent evolution possibly related to a 2-to-l dimensional crossover in the be- 
havior of the system. From a comparison of O-poor vs O-rich PZT films, it appears that 
oxygen vacancies are particularly important for pinning domain walls, and allow even rela- 
tively small domains to be stabilized at high temperatures, close to To This observation is 
especially interesting for potential applications, for which, in fact, it may be more useful to 
have higher oxygen vacancies and thus more stable domain structures. 
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Ar:0 2 = 58:42, 180 mTorr and Ar:0 2 = 80:20, 225 mTorr were used by P. P. and X. H., 
respectively 

(A) The 50 nm film was heated for 3 hour intervals in 100°C steps to 485° C, then in 50° C 
steps till domain disappearance at ^785°C. The 52 and 91 nm films were heated for 1 hour 
intervals in 100°C steps to 485°C, then in 25°C steps till domain disappearance at ^625°C. 

(B) The 50 nm film was heated for 30 minute intervals in 100°C steps to 300°C, 50°C steps to 
450° C, then 25° C steps till domain disappearance at ^525° C. The 40 nm films were heated for 
20 minute intervals in 100°C steps till 400°C, 50°C steps till 550°C, then 25°C steps till domain 
disappearance at ^625°C. The 100 nm film was heated for 20 minute intervals in 100°C steps 
till 200° C, then 25° C steps till domain disappearance at ^325° C. 

The significant asymmetry in domain switching and stability as a function of the polarity of the 
tip bias is a feature already noted in previous studies 44 and can be related to the particular device 
geometry used. Essentially, the metallic electrodes on either side of the ferroelectric thin films 
can be considered as two back-to-back Schottky diodes, with the barrier height determined by 
the work functions of the respective electrode material d 40 * 45 *. In addition, further asymmetry in 
current- volt age and switching characteristics can result from strain and other boundary-specific 
conditions, and, in thin films, the direction of polarization. Finally, the device architecture itself 
is strongly asymmetric: one electrode is planar, while the other is the nanoscale AFM tip, 
with a nominal radius at apex form metal-coated tips given as 25-50 nm. In such devices, the 
nucleation and growth of domains is therefore strongly dependent on the polarity of the voltage 
applied to the AFM tip. 
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